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Abstract 

Abstractions of crystalline materials known as periodic body-and-bar frameworks are made of rigid 
bodies connected by fixed-length bars and subject to the action of a group of translations. In this paper, 
we give a Maxwell-Laman characterization for generic minimally rigid periodic body-and-bar frameworks. 
As a consequence we obtain efficient polynomial time algorithms for their recognition based on matroid 
partition and pebble games. 

1 Introduction 

In this paper, we study dimensional periodic frameworks made of rigid bodies connected with rigid bars. 
Fig. 1 gives an example. We prove a combinatorial characterization for the quotient graphs of generic 
minimally rigid frameworks, in terms of matroid unions of graphs satisfying Maxwell-sparsity conditions. 
As a consequence, we obtain an efficient, polynomial time algorithm for their recognition, based on matroid 
partition and pebble games. 




Figure 1: Fragment of a periodic body-and-bar framework associated to Prussian Blue crystalline materials [11]. 

Crystalline materials. Periodic structures are naturally associated with crystallographic studies [ ]. The 
central role played by a periodicity lattice and corresponding fundamental domains was recognized even before 
X-ray diffraction opened up experimental possibilities. Crystals known as framework materials [ , 26, 27] 
have distinctive substructures which can be modeled and investigated as articulated systems made of rigid 
pieces and joints described abstractly by an infinite periodic graph. Similar models appear in structural 
engineering in connection with infinite trusses, foams or cellular materials [j, 7, 12, 16]. While a fair 
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number of traditional areas of mathematics are related or dedicated to lattice structures and periodicity, 
from crystallographic groups and quadratic forms to the geometry of numbers and sphere packings [5, 58], 
mathematical studies of periodic frameworks are of relatively recent date [ , , ]. 

Displacive phase transitions. Crystallographers have long observed that under variations of temperature 
or pressure, materials may undergo structural changes while retaining their crystalline nature. Understanding 
these phase transitions remains a challenging problem [ ]. The role of geometry is particularly relevant in 
displacive phase transitions which may be interpreted as continuous deformations of a given framework [ ]. A 
rigorous mathematical understanding of rigidity and flexibility properties of periodic frameworks may assist, 
clarify and complement approaches based on computational physics and materials science [ , 18]. 

Computational challenges. Flexibility studies of macromolecules modeled as mechanical (finite or peri- 
odic) frameworks [ , ] are oftentimes conducted with computationally expensive yet numerically imprecise 
simulations. For finite structures characterized by theorems of Maxwell-Laman type (described below), much 
faster approaches based on degree-of-freedom counting and rigid component calculations are known. But such 
theorems are exceedingly rare and difficult to obtain. For periodic structures, one would expect, intuitively, 
that the number and distribution of connections in a "unit cell" would dictate the flexibility properties of the 
framework, but the occurrence of certain 'paradoxes' [18] related to choices of periodicity has stymied the 
progress. Making precise such an intuition with an appropriate mathematical definition has emerged only 
very recently [ ], opening the way to a combinatorial treatment and efficient algorithms [ ] for rigidity and 
ffexibility analysis of periodic bar- and- joint structures (made from rigid bars connected through rot at able 
joints). 

The approach adopted in this paper is rooted in the theory developed in these recent papers [ ', ^], where 
we gave a Maxwell-Laman theorem for generic periodic bar-and-joint rigidity, thus reducing both the long- 
standing ffnite case and several other periodic situations to genericity reffnement conjectures. It is therefore 
important to understand in this context the intrinsic theoretical difficulties we were facing in proving the 
results of this paper: the rest of the introduction gives this perspective. 

Modeling crystalline frameworks. Because of interatomic bond length and angle constraints, some 
crystalline materials may be abstracted as a periodic framework made from rigid bars connected at vertices, 
to which, in principle, the theory developed in [ , ] applies. Often, substructures making individual rigid 
bodies are identified from the outset, such as the octahedra in Fig. 1 or, in a 2D example, the gray bodies 
from Fig. 2(b). In this paper we focus on those made of rigid bodies connected by rigid bars, as illustrated 
in Fig. 1 and 2(b). 




Figure 2: (a) A 2D periodic bar-and-joint framework containing smaller rigid components. In (b), the same 
framework, viewed as a body-and-bar framework: the rigid components form the bodies, and the remaining bars 
connect distinct bodies. 

The practical difficulty in applying the results of [ ] to these frameworks comes from the generic nature of 
this characterization theorem. Indeed, the placement of the "bars" in actual crystal models, with several 
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bars sharing endpoints and connecting pairs of atoms with specified regularity, is not guaranteed, a priori^ 
to be generic. This is a very subtle distinction for which the reader may need further guidance. Hence 
we discuss now briefiy the challenges inherent in addressing questions of generic rigidity, at the same time 
presenting an overview of the techniques we use. 

Maxwell-Laman sparsity conditions. The theory of finite frameworks has a distinguished tradition going 
back almost 150 years to Maxwell [ ], where a sparsity condition was shown to be necessary for minimal 
rigidity of bar-and-joint frameworks: in dimension for any subset of d < < \V\ vertices, the graph 
should span at most dn' — (^^^) edges, with equality for the whole set of n = |V| vertices. Its sufficiency 
for generic frameworks in dimension two was proven over 100 years later (Laman's theorem [ ]), and is 
known to fail in higher dimensions. The problem of completing the combinatorial characterization for bar- 
and-joint frameworks in arbitrary dimensions remains the most conspicuous open question in rigidity theory. 
However, some restricted classes of finite frameworks have been shown to have similar Maxwell-Laman 
counting characterizations: body-and-bar and body- and- hinge frameworks [17, 33, 34, 35, 39], panel-and- 
hinge frameworks [ ] and {d — 2)-dimensional plate-and-bar frameworks [33, 35, 36]. 

Generic structures. All known theorems giving combinatorial characterizations of framework rigidity and 
fiexibility proceed by associating a rigidity matrix M to the framework. A set of algebraic constraints (e.g. 
fixed distances between some pairs of points) are applied on a collection of rigid objects (e.g. a finite set 
of points), with the goal of obtaining a rigid structure, i.e. one in which all the pairwise distances are 
determined. The tangent space at a generic point on the variety given by these algebraic constraints induces 
a matrix, called the rigidity matrix^ whose rank we seek to compute. The maximal number of independent 
constraints that can be imposed is also the minimum number that would make the structure rigid. Together 
with additional constraints that eliminate the trivial rigid motions, they lead to a square rigidity matrix. 
Maximal rank, attained when the determinant det{M) of the rigidity matrix is non-zero, indicates that the 
framework is minimally rigid: it is rigid and the removal of any constraint will make it fiexible. The condition 
of maximal rank is therefore an algebraic condition in the free variables of the system. A framework (with 
a square rigidity matrix) is generic if its actual geometric realization falls outside the set of zeros of the 
algebraic equation det{M) = 0. This extends naturally to frameworks whose rigidity matrix is not square, 
by expressing the rank condition in terms of corresponding minors of M. 

Theorems of Maxwell-Laman type for minimally rigid structures characterize classes of graphs underlying 
frameworks whose rigidity matrices have maximum rank, generically, for certain structural models as de- 
scribed above. A specific geometric realization of a graph belonging to such a class is, with probability 1, 
rigid. But the measure-zero set when it is not rigid is non-empty, and may contain important cases that 
appear in practice. Refining the genericity conditions to prove that they include such practical occurrences 
is in general a very difficult problem. A remarkable and important success of this nature is the recent proof 
of the Molecular Conjecture [ ]. Our contribution here is of a similar fiavor, in the sense that it eliminates 
a certain type of genericity assumption. 

Maxwell-Laman sparsity for periodic frameworks. The connection pattern of a periodic body-and- 
bar system determines an infinite multi-graph G = (V^E)^ with vertices corresponding to bodies and edges 
corresponding to bars between pairs of connected bodies. See Fig. 2(c) for an example in 2D and Fig. 1 for 
one in 3D. Periodicity, or more precisely d-periodicity (where d represents the dimension of the ambient space 
in which the graph is realized geometrically) requires a free Abelian automorphism subgroup F C Aut{G) of 
rank d. We work under the assumption that the quotient graph G/F has a finite number n of vertex orbits 
and a finite number m of edge orbits. 

The problem of characterizing periodic frameworks is dramatically different from the finite case: the generic 
periodic bar-and-joint frameworks have been characterized in all dimensions [ ] by a Maxwell- sparsity con- 
dition on the quotient graph. The finite bar-and-joint frameworks, as well as the periodic body-and-bar, 
body-and-hinge, body-and-pin etc. appear as special cases where additional algebraic dependencies are 
present. As indicated, such cases are not guaranteed to be generic a priori, and even getting a necessary 
sparsity condition (something that was trivial in the finite case) is not easy. 
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Matroidal sparsity conditions. An additional difficulty, relevant from an algorithmic point of view, arises 
from the fact that not all types of sparsity conditions are matroidal. For instance. Maxwell sparsity for finite 
3D bar-and-joint frameworks (the "3n-6" condition) is not. For matroidal conditions, efficient recognition 
algorithms exist. Relevant for our setting are matroid partition [ ] and pebble game algorithms [22]. 

We can now state our contributions in more technical terms. 

Contribution and novelty. Motivated by foundational questions in theoretical crystallography, we present 
a complete mathematical and algorithmic solution to the periodic body- and-bar problem by giving a Maxwell- 
Laman characterization. We formally define the periodicity model, identify the necessary sparsity conditions, 
prove their sufficiency and, finally, give efficient algorithms for deciding rigidity, based on matroid partitioning 
and pebble games. 

As already seen in the bar-and-joint case [1, 2], the mutual relations between the finite and the periodic 
settings are intricate. We emphasize here that our approach considers the periodicity group F as part of 
the initial data (G, F) but allows the variation of its representation as a lattice of translations of full rank. 
Moreover, the problem of minimal rigidity may be formulated at several levels. The most basic level looks 
only at the structure of the quotient graph G/T and provides characterizations 'up to a generic lifting of 
edges' from G/F to a covering periodic graph. The main results presented below are formulated at this level. 
For n = 1, that is, in the case of a single body orbit under F, an answer is given at all levels and offers a 
first measure of the contrast between solving in terms of G/T or fully in terms of the given periodic graph 
(G, F). We discuss aspects implicated at this higher level after solving the problem at the basic level. 

Related work. Among other recent contributions to the subject, we mention [24] which addresses the 
two-dimensional case of periodic bar-and-joint frameworks and [29, 3' ] where certain classes with 'added 
symmetry' are investigated. The final part of the thesis [ ] includes some considerations about body-and- 
bar frameworks "on the fixed torus". This setup corresponds with a fixed lattice of periods. A necessary 
condition for minimal rigidity is presented and conjectured to be also sufficient. It will be seen below that 
this is a restricted version of the minimal rigidity problem in terms of (G,F). 

Overview. We describe our setting in Section 2, where we also derive the rigidity matrix. We indicate 
different levels for the problem of minimal rigidity in Section 3 and solve it at all levels for the base case 
n = 1. Our main result characterizing minimal rigidity in terms of the structure of the quotient graph G/F is 
proven in Section 4. We then extend these results in Section 5 to mixed periodic plate-and-bar frameworks. 
We conclude with the algorithmic consequences and several open problems. 

2 Periodic graphs and frameworks 

In this section we introduce our formal definition of periodic rigidity for periodic body-and-bar frameworks. 
For standard rigidity theoretic definitions such as infinitesimal rigidity and rigidity matrix, see [., 39]. 

We remind the reader a classical concept from finite rigidity theory: a finite body-and-bar (i-framework is 
a finite collection of rigid bodies in connected through bars. The endpoints of a bar lie on two distinct 
bodies and act as rotatable joints, i.e. they permit rigid motions of the incident objects (bodies or bars) 
relative to each other, while contact is maintained at the joints. Several bars may be placed between the 
same pair of bodies. See Fig. 3. The combinatorial (incidence) structure of such a framework is captured 
by a multi-graph G = {V^E)^ whose vertices V correspond to bodies and edges to bars. Body-and-bar 
frameworks have been studied in several foundational papers (e.g. [34, 3^]) and have practical application 
in studies of molecular fiexibility [10, 37]. We proceed now to our new definitions. 

A d-periodic body-and-bar graph is a pair (G, F) of an infinite multi-graph G = (V, E) and a group 
F acting on it. The group F, called the periodicity group of G, is a rank d free Abelian subgroup of the 
automorphism group Aut{G) of G, which acts without fixed points and has a finite number of vertex and 
edge orbits. The elements 7 G F may be called periods of G. Since both n = card{V/T) and m = card{E /T) 
are finite, this setting induces a finite quotient multi-graph G/F = (V/F,£^/F). 
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(a) 

Figure 3: A finite, minimally rigid body- and- bar framework in 2D, with 3 bodies and 6 bars, and its associated 
multi-graph. 

Isometries. We denote by E{d) the Euclidean group in dimension that is, the isometry group of the 
Euclidean space R^. The subgroup of translations is denoted by T{R^) C SE{d). The connected component 
of the identity, denoted by SE{d), is made of all the orientation preserving isometries and is referred to as the 
special Euclidean group or the group of rigid motions in R^. A transformation in E{d), respectively SE{d), 
can be represented by a pair T = (p, M), where p G R^ denotes a translation and M G 0((i), respectively M G 
SO{d)^ is an orthogonal transformation, resp. special orthogonal one. M is a matrix with \det{M)\ = 1, 
resp. det{M) = 1. An Euclidean transformation T : R^ ^ R^ operates by the formula T{x) = Mx -\- p. In 
this notation translations correspond to pairs (p^Id)^ where Id is the dx d identity matrix. 




(a) 

Figure 4: Top: A small fragment of an abstract periodic graph. Only the edges between depicted vertices are drawn. 
Equivalent vertices under the action of the group of periods have the same label. The quotient multi-graph and a 
realization are shown in Fig. 6 and Fig. 5. 

To define a body-and-bar framework, we will specify a placement in R^ of the bodies. We assume that all 
our bodies are oriented as the ambient space R^. Then the placement of a body in R^ is described by a 
pair T = {p^M) G SE{d). In other words, we conceive of the body as described by a (positively oriented) 
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Cartesian frame with origin at p and basis vectors corresponding to the columns of the orthogonal matrix 
M eSO{d). 

Our next goal is to define a periodic body- and- bar framework. A 2D example is given in Fig. 5. 

A presentation in of a d-periodic body-and-bar graph (G, F) is given by an assignment of Cartesian 
frames to the vertices r : V ^ SE{d)^ alongside with an injective representation of the periodicity group F, 
TT : F ^ T{R^). We refer to 7r(F) as the lattice of periods and ask that it has rank d. 

Furthermore, for each edge e = (i, j) G E two (arbitrary but fixed) endpoints are given, indicating where 
the bar is attached on the bodies corresponding to vertices i and j. Let us denote these endpoints with 
= q^{e)^ = q^{e). The coordinates G are given with respect to the Cartesian frame marking the 
body corresponding to i, so that for a frame r{i) = {pi^Mi)^ this end of the bar is at the point of R^ with 
coordinates q'^ = Miq^ -\-Pi- 

All this data must respect the conditions required by periodicity, i.e. when we act on a vertex i by a period 
7 G F, we have r(7i) = 7r(7)r(i) and ^*(e) = q^'^{je) etc. The main definition can now be given: 




(a) (b) (c) 

Figure 5: (a) Fragment of a 2-periodic body-and-bar framework (d = 2), with the generators (which are not 
bars) of the periodicity lattice marked by arrows and the orbit of one bar endpoint visualized as a "lattice", (b) 
Representatives of the n = 2 equivalence classes of bodies and m = 6 classes of bars are shown in darker shades. The 
endpoints of the bars, marked in red, are not all distinct in this case, (c) The framework is periodically flexible, as 
illustrated by this deformation of (b). Notice the deformation of both the framework and of the periodicity lattice. 

A d-periodic body-and-bar framework is a connected (i-periodic graph (G, F) together with a presenta- 
tion (r, 7r,g') where all bars have non-zero lengths. 

Realizations, configurations, rigid and flexible frameworks. Presentations which use the same at- 
tachment points q (relative to the corresponding bodies) and induce the same system i of bar lengths form 
the realization space for the data {G^T^q^i). 

Realizations that differ by an isometry of R^ will be considered as the same configuration, and the config- 
uration space of (G, F, r, tt, q) is defined as the quotient space of all realizations by the group E{d) of all 
isometrics of R^. The deformation space of a periodic framework (G, F,r, 7r,g) is the connected component 
of the corresponding configuration. When the deformation space consists of a single point, the framework is 
called rigid, otherwise it is flexible. An example of a fiexible periodic framework is illustrated in Fig. 5, with 
its quotient graph shown in Fig. 6. 

As in [1], this deformation theory has the following useful characteristics: (i) The periodicity group F is part 
of the structure and the deformations of a framework are those preserving this specified periodicity, (ii) The 
lattice of translations 7r(F) representing F is allowed to vary as the framework deforms, (iii) The realization 
space is the solution space of a finite algebraic system of quadratic equations, (iv) There is an infinitesimal 
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(a) 

Figure 6: Top: representatives of orbits for bodies and bars for the periodic framework from Fig. 5. Bottom: the 
quotient multi-graph. 

deformation theory, obtained by differentiating the equations. This leads to the concept of infinitesimally 
rigid framework^ characterized by a rigidity matrix of maximum rank. 

In particular, standard arguments from algebraic geometry (inverse function theorem) can now be used to 
prove that infinitesimal rigidity implies rigidity. 

Minimal rigidity. Following a heuristic going back to Maxwell [ ], we do a quick calculation of the 
minimal number of parameters needed to specify a rigid framework. This gives us an indication of the 
number of bars we would expect in a minimally rigid graph. A vertex orbit requires {^~^^) parameters to be 
specified, and an additional d'^ parameters specify the periodicity lattice. The total number of variables in 
the algebraic system is thus ^(^^^) + d'^. Each edge orbit leads to a constraint (equation), and generically 
it eliminates one degree of freedom. Since we cannot eliminate the trivial isometrics, we expect to have at 
most -\- d'^ — {^~^^) = (n — 1)(^^"'^) + d^ independent equations. By definition: 

A minimally rigid d-periodic body-and-bar graph (G, F) is a d-periodic graph with m = (n — l) (^^^) + 
d^ which has an infinitesimally rigid presentation. 

For n = 2 bodies in dimension 2 this number is 7. Since the quotient graph from Fig. 6 of the periodic 
framework from Fig. 5 has 6 < 7 bars, we expect the periodic framework to be flexible with one degree of 
freedom (Idof). 

We have now all the ingredients to start developing the rigidity theory of periodic body-and-bar frameworks. 
The first step is to define and analyze the rigidity matrix associated to a generic framework. 

The rigidity matrix. We express the rigidity matrix in terms of the following coordinates: we choose a 
basis for the lattice of periods 7r(r): Ai,...,Ac^ and denote by A the dx d matrix with columns A^; then 
we choose representatives for the n equivalence classes of vertices modulo F and consider the corresponding 
bodies as marked by {pi^Mi)^ i = l,...,n. We then choose an orientation of the edge set in the quotient 
graph G/F and select representatives for the equivalence classes of edges modulo F by the rule that all 
edge representatives originate at the already chosen vertex representatives. We remind the reader that the 
endpoints q'^^ of an edge representative are fixed parameters, but their choice is only subject to q^ q^ - In 
order to keep the notation as simple as possible, we do not pursue the details related to marking distinctions 
between several bars possibly joining the same pair of bodies. 

Since one end (the origin) of any edge representative is on a body representative, the other end is on the 
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corresponding representative translated by a period A = A(7(e) = Ac, where c = (7(e) is a column vector 
with integer entries i = 1, (i. 

Then, a bar constraint corresponding to an edge e between bodies i and j, with endpoints g'*(e), resp. q^{e), 
on these two bodies, takes the form ||AC(e) + Mjq^{e) — Miq^{e) -\- pj — = i^, with e running over all 
m = \E/T\ edge representatives. 

We will eliminate the equivalence under rigid motions by assuming that the first representative body is fixed, 
with (pi,Mi) = (0,/(i). We note here that for an orthogonal matrix M, an 'infinitesimal variation' M takes 
the form M = MA^ with A skew-symmetric. Thus, if we differentiate the system at (p2, ^2), (Pn, ^n), A, 
we obtain a linear system in infinitesimal variations (p^, ^4^), i = 2, n and A with rows 

(Ac + MjAjQ^ - MiAiq' + pj - Pi, Ac + MjQ^ - Miq' + pj - Pi) = (1) 

where we have dropped the notational dependency on e and used indices i, j for the chosen body-representatives, 
with Ml = Id-, Pi =0 and Ai = 0. The rigidity matrix is the matrix of this linear system (for some specific 
ordering of the unknowns). 

Normalized coordinates. Given any (i-periodic body-and-bar framework, the manner of marking chosen 
representatives in each vertex orbit by Cartesian frames (p^. Mi) is also a matter of choice and we may adopt 
the normalization which has all = and Mi = Id- The rigidity matrix is then the matrix of the linear 
system in pi^ Ai = —A\ and A: 

(Ac + AjQ^ - Aiq' + Pj - Ac + - q') = (2) 
where the edge vectors Ac + q^ — q^ run over a complete set of edge representatives. 

3 Minimal rigidity 

As already observed, up to Euclidean rigid motions, a periodic body-and-bar framework is described by 
— 1) + (i^ parameters. A bar equivalence class imposes a single constraint on these parameters 
hence rigidity requires at least — 1) + constraints. A (i-periodic body-and-bar graph (G,r), with 

n = \E/T\ and m = \E/r\ = (^^^)(^ — 1) + d^, is called minimally rigid if it can be presented as an 
infinitesimally rigid framework in R^^ that is, as a framework with rigidity matrix of maximal rank m. 

Liftings and marked quotient graphs. The problem of characterizing the structure of minimally rigid 
periodic graphs can be formulated at three distinct levels^ illustrated in Fig. 7. The most discerning and 
demanding level concerns the periodic graph (G, F). Two other less discerning levels focus only on the 
quotient multi-graph G/F, with or without some additional lifting information retained. 

The information lost upon passing from the periodic framework (G, F) to the quotient graph G/T is the lifting 
of the edges to specific vertex representatives of the same orbit. This information can be retained in the 
following form. We fix a basis of F, that is, an isomorphism T ^ Z^. Then we choose vertex representatives 
and an orientation of the quotient graph. By using edge representatives which originate at the chosen vertex 
representatives, we have a well defined lifting function C : E/T ^ which allows retrieval of (G, F) from 
G/F. When several bars connect the same bodies, we can retain this information by joining the vertex 
representatives with single edge, marked with a number indicating the edge multiplicity, as in Fig. 7(b). 
This 'multi-graph with multiplicities' is denoted by Q/T. Combining the edge multiplicities with the lifting 
function yields a labeled multi-graph with multiplicities, such as illustrated in Fig. 7(c). This contains all 
the information needed to reconstruct (G,F). 

Remark: An invariant way for describing the lifting data can be given in terms of the first homology group 
(with integer coefficients) i^i(G/F), where G/F is now seen as the topological space corresponding to the 
one-dimensional C W-complex defined by the quotient graph [ ] . With G considered in the same topological 
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(a) (b) (c) 



Figure 7: (a) The basic quotient graph G/T. (b) The quotient graph with multiphcities Q/T. (c) The 
labeled and oriented quotient multi-graph with multiplicities, equivalent with (G, F). 

perspective, we have a (connected) Abelian covering G G/T with covering group F, and this is equivalent 
to a surjective group homomorphism Hi{G /V) F. 

In the rest of this paper, we give a complete characterization of (G, F) for the case of a single body n = 1, 
and a characterization in terms of the quotient graph G/T for arbitrary n. 

The n = 1 case. The characterization of minimal rigidity in terms of G/T is trivial for n = 1, since the 
only condition about the quotient graph is that of having cP loops. This is enough for finding generic liftings 
which have infinitesimally rigid framework presentations. We address here the more refined questions related 
to (G,F) and ^/F, assuming m = d? . 

The quotient multi-graph with multiplicities Q /T has a single vertex and a number of loops with multiplicities 
ki^ z = 1, with Yl\=i = d^- 

Minimal rigidity in this setting is equivalent with the condition that all multiplicities ki be at most d. This 
will follow from the solution of the sharper problem about the structure of (G, F). 

Recalling that n = 1, m = and that we eliminate the equivalence under rigid motions by fixing (pi. Mi) = 
(0, Id), we obtain from (1) a rigidity matrix with d'^ rows of the form G(e) [AC{e) + {q^ {e) — g*(e))], where 
the tensor notation stand for a listing of all the d'^ products of the d components of the two indicated vectors. 
Since the endpoints of bars can be chosen at will, the edge vectors AG(e) + (<7^(e) — <7^(e)) are arbitrary 
(non-zero) vectors, say h = h{e). With simplified notation q for the coordinates in G(e), such a row can be 
presented as [ hi{ci C2 ... q) /i2(ci C2 ... q) ... hd{ci C2 ... q) ]. 

In this form one recognize a union operation on d copies of a linear matroid, leading to 

Theorem 1. Let (G, F) be a d-periodic body-and-bar graph with n = \V/T\ = 1 and m = \E/T\ = d^ . Then 
(G, F) has infinitesimally rigid framework presentations in if and only if for any set of equivalence classes 
of edges F C E/T we have 



\F\ <d-dim{spanC{F)) (3) 

Proof: The necessity follows directly from the form, described above, of the row constraints. The sufficiency 
follows from results in matroid theory, which will be used again in arguments for subsequent statements. 
Specifically, d - dim{spanC (F)) is an increasing submodular function with values in the non-negative integers 
[28]. By [30], the associated matroid is the union of d copies of the matroid associated to the function 
dim{spanC {F)) . Hence a basis of the former is a disjoint union of bases in the latter which yield the 
maximal rank of the rigidity matrix. □ 

This proof also implies our earlier claim that a quotient graph Q /T comes from a minimally rigid (S,F) if 
and only if all multiplicities ki, i = 1, are at most d. 
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Indeed, any choice of lifting function C : E /V Z^, with image consisting of classes of at most d vectors 
(corresponding to the given multiplicities) and with representatives in general position, will satisfy the 
criterion of Theorem 1. 

4 Minimal rigidity characterization of the quotient graph G/T 

In this section we present the main theoretical result. We assume that m = (^^^) {n — l)-\-(P and characterize 
the quotient graphs G/T which are obtained from minimally rigid (i-periodic body-and-bar graphs (G, F). 
In other words, we characterize minimal rigidity 'up to a generic lifting of the edge representatives'. 

Graph sparsity. We say that a multi-graph G = (V^E) (possibly with loops) is (a, 6)-sparse, or has an — b 
sparsity type, if any subset of < |V| vertices spans at most an' — b edges. When equality holds we have a 
tight (sub)graph. See [ ] for a comprehensive treatment of this kind of graph sparsity. 

When < b < 2a, in particular when 6 = and b = a, the family of (a, 6)-sparse graphs on a fixed number of 
vertices n form a matroid. When b = a, this matroid is the union of b graphic (tree) matroids. When 6 = 0, 
it is the union of a bi-cycle matroids, whose bases are the (l,0)-sparse graphs (also known in the literature 
as map-graphs, pseudo- forests or bi-cycles). We also obtain a matroid when b < 0: as shown in Corollary 4 
below, it is the union of the (a, 0)-sparsity matroid and the uniform matroid on b elements, or, equivalently, 
it is any graph which has an -\- b edges and contains a (a, 0)-sparse subgraph. 

We characterize now the quotient graphs of periodic minimally rigid graphs in terms of matroid unions of 
graph sparsity matroids. 

Theorem 2. A quotient graph G/T with m = — 1) + is the quotient graph of an infinitesimally 

rigid periodic body-and-bar framework in if and only if it satisfies one and hence both of the following 
equivalent conditions: 

(i) it decomposes into a disjoint union of two graphs on the full set of n vertices, one with dn — d edges and 
sparsity type, the other with (^n + (^X^) edges and sparsity type; 

(a) it contains the disjoint union of two graphs on the full set of n vertices, one with dn — d edges and 
sparsity type, the other with {^n edges and sparsity type. 

As shown in ["] , graphs with dn — d edges on n vertices and with dn — d sparsity type are minimally rigid 
graphs for the problem of rigidly connecting with bars a set of n translated bodies (i.e. parallel Cartesian 
frames). Thus, the presence of a dn — d sparse subgraph is related to the fact that vertex-orbits (when 'frozen 
solid' together into a body) provide a system of n translated bodies which has to become rigid. 

Proof of sufficiency: We begin with some combinatorial background about sparsity of type an + 6, with a 
and b non-negative integers. Again, in this discussion, graphs may have multiple edges and loops. 

Loop breaking. When a graph has loops, we shall describe as loop breaking the operation of replacing a 
loop by an edge connecting the loop vertex with some other vertex of the graph. It is immediate that loop 
breaking preserves an -\- b sparsity. The following lemma is the relevant converse. 

Lemma 3. Let G be a graph with an-\-b edges on n vertices, with a and b non-negative integers. Suppose G 
is an-\-b sparse, that is, any subset of n' vertices has at most an' + b edges between them. Then, there is a 
graph G, made of an-\-b loops on n vertices, which is an-\-b sparse and yields G after an adequate sequence 
of loop breaking. 

Proof: We claim that we can always replace an edge of G between vertices u and v by a loop at one of these 
two vertices so that an -\- b sparsity is preserved. Suppose this were not the case. Then there are subsets of 
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vertices U and with u in U but not in V and v in V but not in U with a\U\ -\-b edges in U and a\V\ -\- b 
edges in V. Since the union U UV has at most a\U UV\ -\- b edges, the intersection U OV must be non-void 
and tight. Even so, the union U UV has aheady a\U UV\ -\- b edges without counting the edge we started 
with. This contradiction proves the claim and the lemma follows by iteration. Note that the resulting graph 
G must have at least a loops at each vertex. □ 

Corollary 4. With the notations of the above lemma, G is an-\-b sparse if and only if it contains a subgraph 
with \V\ = an edges on the n vertices and sparsity type an. 



We show now that a loop breaking operation may be performed on minimal infinitesimally rigid frameworks. 
Let us assume that such a framework has a loop at vertex j in the quotient graph G/T. Using normalized 
coordinates as in (2), the loop corresponds to an equation of the form: 

(Ac + Aj{q^ - q'^), Ac + - q'^) = (4) 

which involves only the difference q^ — q'^ . This allows the assumption q'^ = for the origin of the bar. We 
are going to replace this loop at vertex j of G/T by an edge from vertex i to vertex j. The corresponding 
bar has one endpoint at the origin of the body representative indexed by i and the other endpoint at kq^ on 
the body representative indexed by j translated by period A/cc, for some integer k. This gives an equation 
of the form: 

{Akc + Ajkq^ + pj - Pi, Akc + kq^) = (5) 
which, divided by /c^, takes the form: 

(Ac + Ajq^ + ^{pj - Pi), Ac + q^) = (6) 

With all other rows unchanged, the rigidity matrices for the two frameworks involve only the exchange of row 

(4), where q'j = 0, with row (6) and for a sufficiently large integer /c, the rank won't be affected. This shows 
that loop breaking operations can be performed on quotients of minimally rigid graphs with preservation of 
minimal rigidity. 



The sufficiency proof follows from these results. Indeed, if G/T has edges separated into di dn — d tight sparse 
graph and a (2)^ + (^2^^) tight sparse graph, we can produce minimally rigid d-periodic graphs with that 
quotient by a two-step procedure. First, we use the dn — d sparse subgraph, but replace the remainder graph 
with a graph made only of loops, as guaranteed by Lemma 3. Infinitesimally rigid frameworks corresponding 
to this quotient structure can be constructed in a natural way, as indicated immediately below. Then, as a 
second step, we break the appropriate loops and obtain the initial graph while maintaining minimal rigidity. 

When the (2)^+ C^2^) tight sparse subgraph is made entirely of loops, natural liftings of G/T to minimally 
rigid (G, F) are constructed as follows. First, we use {^~^^) loops to make rigid the lattice of periods, as in 
[1, 2]. Note that bars placed as periods can be interpreted as loops at any vertex of the quotient graph. 
We distribute them so that exactly (2) loops remain for bar assignment at each vertex class. Now in each 
vertex orbit we place the available (2) classes of bars so that the whole body-and-bar orbit becomes a single 
rigid body. Thus, with the (2)^+ (^t^) ^^^P subgraph we have turned the system into n parallel bodies and 
now the dn — d sparse subgraph is precisely what we need to turn this system into a single rigid body, as in 
[2]. □ 

An alternative proof of sufficiency is given below after completing the necessity. 

Proof of necessity: We consider the m x u rigidity matrix of an infinitesimally rigid periodic body-and-bar 
framework in R^. The number of rows m = (^^^)(^ — 1) + is the number of edge orbits and the number 
of columns u = -\- d'^ = m + (^2^^) is the number of ^unknowns ^ in the linear constraint system for 

infinitesimal variations. By assumption, the rows are linearly independent and the (^^^) -dimensional kernel 
of our rigidity matrix consists precisely of the trivial infinitesimal rigid motions. 
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When we consider some subset of rows, denoted say by F, we have \F\ = rk{F) = u — dim{Ker{F)). 

We denote by the number of vertex orbits incident to the given edge orbits and hy Uf the number of 
connected components of F as a set of edges in the quotient graph. We observe that 

dim{Ker{F)) >(^^^^'^+{n- np) (^^^^^+ dicop - 1) (7) 

where the numbers on the right hand side count the fohowing independent infinitesimal deformations re- 
specting the constraint system F: trivial infinitesimal rigid motions, arbitrary infinitesimal deformations 
for the body representatives for non-incident vertices and relative infinitesimal translations of the connected 
components in F. Thus, we have a necessary edge sparsity condition: 

\F\=u- dim{Ker{F)) < d{nF - ujf) + Q + ("^ J ^) 



Written in this form, the right hand side becomes very suggestive from the point of view of matroid theory 
[28]. Indeed, if we consider a ground set of sufficiently many edges and loops on n vertices (so that any 
graph with m edges can be found as a subgraph), then the function defined by the right hand side is an 
increasing submodular function with values in the non-negative integers and determines a matroid which 
is, by a theorem of Pym and Perfect [ ], the union of d copies of the graphic matroid determined by 
fi{F) = riF — ojf-, then copies of the hi-cycle matroid determined by /2(F) = uf and the uniform 
matroid determined by /3(F) = (^^"^)- 

A basis in a union of d graphic matroids is a union of d spanning trees and this is precisely what dn — d tight 
sparse graphs are, while (2)^ + (^2^^) tight sparse graphs give the bases of the remaining union. Thus, the 
m edges corresponding to our maximal rank rigidity matrix have the claimed decomposition. □ 

This completes the proof of our main Theorem 2. 

The result proven above for minimally rigid d-periodic body-and-bar frameworks is obviously related to 
the characterization obtained in [2] for minimally rigid d-periodic bar-and-joint graphs. These two types of 
periodic framework structures may be seen as the extreme cases k = d^ respectively /c = 0, of the family 
obtained by articulating with bars /c-dimensional plates, with < k < d. In the full paper, we show that our 
proof technique can be extended to these general periodic frameworks. 

To illustrate other techniques, we present now an alternative proof of sufficiency for Theorem 2. 

Suppose the bar corresponding to e G E connects two bodies Bi and Bj^ which are identified with Cartesian 
frames {pi^Mi) and (pj^Mj) with Pi^Pj G and Mi^Mj G SO{d). We normalize them (i.e., = and 
Mi = Id). As we explained previously, a bar gives a linear equation 

(Ac + AjQj - AiQi + pj - Ac + Qj - Qi) = 0, (9) 

where Ai is a skew-symmetric matrix with Mi = MiAi and Qi^qj are endpoints of the bar in each frame. 

Let Wi be a vector in which is obtained by aligning (2) independent entries of Ai. By setting h = 
Qj -\- Ac — Qi, we can simply write (9) by 

{pj -Pi,h) + {{wj.qj Ah) - (wi.qiAh)) + (Ac, /i) = 0. (10) 

The rigidity matrix is then redefine as the matrix of the linear system (10) in Pi^Wi and A over all edge 
representatives. The columns of the rigidity matrix are naturally divided into three parts, the translation- 
part, the rotation-part, and the lattice-part. Namely, they are the column spaces associated with w and 
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A, respectively, and have the dimensions (in, (2)^5 {^~2^)^ respectively. Also each row of the rigidity 
matrix is a {dn + (2)^ + d^) -dimensional vector whose coordinates can be also divided into three parts. 

Let be the i-th element of the standard basis of R^. The following proposition gives an alternative proof 
of the sufficiency of our main theorem. 

Proposition 5. Suppose the set of edge representatives can be partitioned into three parts: 

• d edge- disjoint spanning trees Ti^. . .Td, 

• (2) edge-disjoint spanning pseudo-forests Fij for 1 < i < k and i < j < d (where we may assume that 
in each Fij the edge set is oriented so that every vertex has in- degree exactly one), and 

• U = {e1.:l<i<d,l<j <i}. 
Define a realization by 

• ^ = h, 

• C(e) = e^, qj = = for e e Ti, 

• C(e) = ej.qj = e^, qi = for e e Fij, 

• C{e*j) = N{ei + ej),qi = qj = for e*^ with 1 < i < d and i < j <d, 

where N is an integer with N > n. Then, the rigidity matrix takes full rank. 

Proof: Recall that A = denotes the infinitesimal motion of axes of the lattice. Since any framework has 
(^^^) trivial motions, we shall eliminate such trivial motions by restricting translations of a body by pinning 
down a point and by restricting rotations of the lattice. Namely, let py* = for a vertex v"" and \ij = for 
every (i, j) with 1 < i < d — 1 and 1 < j < i. To do that, we add new (^J^) rows to the rigidity matrix in 
an appropriate manner. We now prove the resulting matrix has full rank. 

We first note the entries of row for each e G Fij. Since ^2 A /i = ^2 A (^2 + C'(e) — ^1) = A and gi A /i = 0, 
the row in the rotation-part is written as 

(0, --^O, 1 ,0, --^O) 

for each e = uv ^ Fij (where we indexed each coordinate of a vector in by a triple {v' ^i^j) with v' 

and 1 < i < j < d). Thus the set of row vectors of Ui<i<j<d-^u forms the identity matrix of size i^n in 
the rotation-part of the rigidity matrix. Moreover, since ^1 = ^2 = for all edges e in (|J^ Ti) U t/, all entries 
of the rest of rows are zero in the rotation part. 

We thus need to show that the rest of rows vectors are linearly independent when restricting to the column 
space of the translation-lattice-part. We again note the entries of each row: For each e = uv e Ti^ the row 
in the translation part is 

(n,z) ... {v,i) 

(0, --^O, -1,0, •••,0, 1 ,0, --^O), 
while in the lattice part the row has the form (8) e^, that is. 

All ••• Aii ••• Xdd 

( , •••,0, 1 ,0, ). 

For e*j = uv E U the corresponding row in the translation part is 

(0, • • •, 0, - AT, • • • 0, - AT, 0, • • •, 0, AT , • • • 0, AT , 0, • • •, 0) 
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while in the lattice part the row has the form A^e^ (g) A^e^, that is, 

All ••• ^i'i' '^^J ^^'^ ••• Add 

( , •••,0,7V^0••• 0,A^^0••• 0,7V^0••• 0,7V^0, ). 

Here is an important observation: Suppose the endpoints of e*^- are u and v in the quotient graph. Then, 
Ti has the path from u to which consists of a sequence of edges ei, 62, . . . , e^. If we sum up the rows 
corresponding to these edges, the translation-lattice-part becomes 

(n,i) ... {v,i) ... All ••• ^ii ••• Add 

(0,---,0, -l,0,---,0, 1 ,0,---,00,---,0, t,0,---, 0), 

where t is the length of the path, which is less than n. Even if we further multiply it by we have 

(n,i) ... {v,i) ... All ••• ^ii ••• Add 

(0, •••,0, -iV,0, •••,0, 7V,0, •••,00,---,0,tiV,0, 0). 

Since N"^ > tN ^ we can safely eliminate all the non-zero entries of the translation-part in the row of e*^ G U . 
It is easy to observe that the submatrix induced by T^, . . . , T^/ together with d additional rows for eliminating 
trivial translations has the rank dn in the translation-part while the submatrix induced by e*^- together with 

(2) rows for eliminating trivial rotations has the rank d^ in the lattice-part. This shows that the rigidity 
matrix has full rank. 

5 Mixed plate-and-bar periodic structures 

The result proven above for minimally rigid d-periodic body-and-bar frameworks is obviously related to the 
characterization obtained in [ ] for minimally rigid d-periodic bar-and-joint graphs. Actually, these two types 
of periodic framework structures may be seen as the extreme cases k = d, respectively A: = 0, of the family 
obtained by articulating with bars /c-dimensional plates in R^, with < k < d. A plate is understood in this 
context as a rigid object marked by a /c-frame in R^. 

Remark: The distinction between articulations made of {d — l)-dimensional plates, also called panels^ and 
those made of bodies, which are d-dimensional plates, is that bars have to be attached at points of the panel. 
Otherwise, any {d — l)-frame has a unique extension to a d-frame with standard orientation. 

We observe that, since /c-frames in R^ can be parametrized by SE{d)/SO{d — /c), the periodic setting with 
variable period lattice will require for minimal rigidity the relation 

m = [d{k + 1) - (^^ + + _ (^rf + 1^ = -d) + idk - (^^ + 1) )n + + 

and the extended version of the theorem is that a quotient graph comes from an minimally rigid d-periodic 
plate-and-bar framework (with all plates of dimension k) if and only if it contains the disjoint union of two 
graphs on the full set of n vertices, one with dn — d edges and sparsity type, the other with {dk — (^^^))^ 
edges and sparsity type. Note that, as in Theorem 2 (ii), there is no condition on the remaining (^^^) edges 
of G/T. 

We may consider a more general setting by mixing various types of plates. The associated d-periodic graph 
(G, r) will have vertices labeled according to the dimension of the plate they represent. This leads to a 
quotient graph G/T with labels, or weights ky G {0, 1, . . . , for each v G V/T. With vertex orbits labeled 
from 1 to n, we write ki, i = 1, ...,n for the corresponding (dimensional) weights. Then, the same type of 
reasoning will establish the following general characterization of minimal rigidity: 

Theorem 6. A vertex weighted quotient graph G/T with 
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is the quotient graph of an infinitesimally rigid mixed plate- and-bar periodic framework in if and only if 
satisfies the sparsity condition: 

\F\ < din, -co,)+Yl [dK - (^%+ (11) 

for every nonempty F C E/T, where Vp is the set of vertices incident to F. 



Remark: The edge sparsity formula (11) replaces formula (8) of the 'pure' top case ky = d. Its role and 
combinatorial unfolding can be treated similarly. Let k'^ = [dky — (^''2^^)]. The function defined by the right 
hand side of (11) determines a matroid which is the union of d copies of the graphic matroid, the uniform 
matroid with rank (^^^), and the matroid determined by /4(F) = X^vgVf ^v- 

The 'loop breaking scenario ' used in Section 3 remains valid for the mixed case and shows in particular that 
any basis of the matroid determined by f^ is obtained from the graph with /c^ loops at each vertex of V/T 
by loop breaking. 

In this sense, 'archetypal' infinitesimally rigid frameworks are obtained by fixing the periodicity lattice 
through (^^^) bars which are periods, then the k'^ loops at the vertex orbit of v are used to eliminate the 
relative motion of the /c^-frame with respect to the lattice and turn the whole orbit of v into a rigid body. 
Finally, the n F-orbits are a sysytem of translated bodies which is rendered rigid by using the remaining d 
spanning trees. An alternative, more detailed argument for the proof of Theorem 6 is the following. 

Proof: The necessity follows from the same argument as Theorem 2 by counting the number of independent 
motions. We show that the sufficiency also follows from the loop breaking argument given in Theorem 2. 

We first provide the corresponding combinatorial part. Let k'^ = [dky — (^''2^^)] for simplicity. By the theorem 
of Pym and Perfect on increasing submodular functions again, the function defined by the right hand side 
of (11) determines the matroid, which is the union of d copies of the graphic matroid, the uniform matroid 
with rank (^^^), and the matroid determined by f^iF) = X^vgVf which is also an increasing submodular 
function with values in non-negative integers. 

To further decompose the last matroid, let Vi = {v e V : ky > £} for £ = 0^ . . . , d and define g£{F) = iVinVrl' 
Observe f^ = ^iZl{d — £)g£, and hence the matroid determined by f^ can be further decomposed into (2) 
matroids. This gives us a decomposition of J^^gv edges into (2) subsets F^i for 1 < £ < d — 1 and 
1 < k < d — £ such that F/^^ is a pseudo- forest spanning only vertices of Vi (i.e., F/^^ is an edge set on Vi 
with \Ff^£\ = \Vi\ such that each connected component contains exactly one cycle). Each F^^i can be oriented 
so that each vertex in Vi has in-degree exactly one, and this orientation implies the following lemma, which 
supplies the converse direction of loop breaking operations. 

Lemma 7. Let G be a graph with "^y^y k'^ edges on n vertices, with k'^ non-negative integers. Suppose 
\F\ < ^y^Y^ k^y for F C E. Then, there is a graph G, made of k'^ loops for each vertex v eV, which yields 
G after an adequate sequence of loop breakings. 

(We should remark that Lemma 7 actually follows from classical Hakimi's orientation theorem [ ] of undi- 
rected graphs directly.) 

We thus convert G/F to a graph which consists oi dn — d tight sparse subgraph, {^^^) edges, and k'^ loops for 
each G As before (^^^) is further converted to (^^^) loops, and we have the base graph. (^^^) loops are 
used to make rigid the lattice of periods. In each vertex orbit we place k'^ classes of bars, corresponding to 
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k'^ loops attached to so that the whole /c^ -dimensional plate-and-bar orbit becomes a single body. Thus, 
using the Xl^^y k'^ + (^2^^) loops we have turned the system into n parallel bodies, and the remaining dn — d 
sparse subgraph makes this system into a single rigid body [ ]. 

The proof is complete by performing a sequence of loop breaking operations we have shown in Theorem 2. 
From above construction, we may assume that each /c^-dimensional plate in the realization is the one deter- 
mined by the first ky column vectors of Id- Then the rigidity matrix is written in the same form: 

(Ac + AjQj - AiQi + pj - , Ac + Qj - Qi) =0 (12) 

over all edge representatives, where Ai is a skew-symmetric matrix with /c^ non-zero entries. Thus we can 
apply the exactly same argument as the proof of Theorem 2 to show that a loop breaking operation preserves 
infinitesimal rigidity. 

6 On minimal rigidity in terms of (G, F) 

The case n = 1 treated in Section 3 illustrates already the contrast between an answer given in terms of 
G/r, which allows a generic lifting of edges to a d-periodic covering graph, and an answer given in terms of 
the periodic graph (G, F) itself. Theorem 1. We discuss here several aspects which must be implicated when 
minimal rigidity is addressed directly on (G, F). We assume G connected and m = (^^^)(^ — 1) + d^. 

We reviewed in Section 3 the fact that the covering G ^ G/T with covering group F can be equivalently 
described through G/F and a surjective group homomorphism i^i(G/F) F. For a given subset F of edges 
in G/F, we shall consider the following numbers or invariants defined in terms of the quotient graph and 
the surjective homomorphism describing its covering (G, F): nF^ooF^dp^dp' ^ where F' denotes a connected 
component of F when considered as a topological subspace of G/F. 

uf is the number of vertices in G/F which are incident to edges in F; 

cjf is the number of connected components of F as a topological subspace F C G/F; 

For each connected component F' C F C G/T we consider the induced maps 



H, (F') ^ Fi (F) ^ Hi (G/r) ^ r (13) 

and define as the rank of the image of Hi{F) in F, while dp' will be the rank of the image of Hi{F') in 
F. Obviously dp' < dp- 

Now, we may formulate a more refined estimate than (7) for the dimension of the kernel of the matrix with 
row constraints determined by F: 

dim{Ker{F)) > (^^^ + + („ - „^) ("^ J ^) + d{uF - 1) + 

where the new terms on the right hand side take account of the fact that when liftings over a connected 
component F' involve only a subspace of the lattice of periods, there are independent infinitesimal rotations 
in the orthogonal complement which are not constrained by F' ^ while the last term refers to the unconstrained 
part of the lattice of periods. 

With |F| denoting the number of edges in F, we have: 
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The question whether this condition is also sufficient for a combinatorial characterization of minimal rigidity 
is left as an open problem. We note that for n = 1, condition (15) becomes (3) since F is connected and 
dp = dim{spanC (F)) . 

7 Algorithms 

We present now some algorithmic consequences of the main theoretical result, Theorem 2. Since our Maxwell- 
Laman characterization decomposes the matroids associated with minimal rigidity of periodic frameworks 
into smaller pieces of well-known (graphic and bi-cycle) matroids, we can simply apply a matroid partitioning 
algorithm (e.g. [ ]) for checking whether a given quotient graph G/T is realizable as a minimally rigid periodic 
framework. 

We can however give a better algorithm, using condition (ii) Theorem 2. It follows from a theorem of [13] 
that a graph is the disjoint union of a — d tight subgraph and a (2)^ tight subgraph if and only if it is a 
(^^^)n — d tight graph. This type of sparsity is matroidal (i.e., the family of tight subgraphs are the bases 
of a matroid on the edges of the complete graph), allowing for the use of the pebble game algorithms of [22]. 
Thus condition (ii) can be checked in O(n^) time. The same algorithm computes the degrees of freedom of 
a generic realization of a given quotient graph, where m may not be exactly (^^^)(^ — 1) -\- d'^. This is done 
by ffist computing a maximum — d-spaise subgraph and then adding at most (^^^) unused edges. 

The number of edges of the resulting graph gives the maximum size of independent equations, its difference 
from — 1) + being the degree of freedom. 

Conclusions 

Motivated by computational studies of crystalline materials, we introduced periodic body-and-bar frame- 
works. We gave a combinatorial characterization of those which are generically minimally rigid, obtaining 
the ffist polynomial time algorithms for recognizing them, and for computing flexibility parameters (degrees 
of freedom) of those which are not rigid. In the full version of the paper, attached as an appendix, we have 
also characterized mixed plate-and-bar periodic structures, and included additional comments on minimal 
rigidity in terms of (G^T). 

Remaining open questions include: (a) characterizing body- and- hinge and (b) body-and-pin periodic frame- 
works, and (c) verifying whether the Molecular conjecture [ ] holds as well in the periodic case. Besides 
intrinsic theoretical interest, such results would have immediate applications, as these structures appear in 
modeling families of natural crystals. 
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